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Abstract: We find the trapped surface for a collision of two sourceless shock waves in AdS 5 
and conclude that such collisions always lead to a creation of a black hole in the bulk. Due to 



holographic correspondence, in the boundary gauge theory this result proves that a thermalized 
medium (quark-gluon plasma) is produced in heavy ion collisions at strong coupling (albeit in M = 4 
super- Yang-Mills theory). We present new evidence supporting the analytic estimate for the time 
of thermalization that exists in the literature and find that thermalization time is parametrically 
much shorter than the time of shock wave stopping, indicating that our result may be relevant for 



X 

description of heavy ion collision experiments. 
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1. Introduction 

The problem of understanding the physics behind thermalization of the medium produced in ul- 
trarelativistic heavy ion collisions is one of the main open questions in heavy ion theory. It has 
become especially important in recent years after hydrodynamic simulations indicated that a very 
short thermalization time of the order of 1 fm/c is required to describe RHIC data [1,2]. Lately the 
problem of thermalization has been studied in the strong coupling framework of the Anti-de Sitter 
space/conformal field theory (AdS/CFT) correspondence [3-5] with the goal of learning about the 
dynamics of the strongly-coupled QCD medium by studying the strongly coupled medium in Af = 4 
super- Yang-Mills (SYM) theory [6-14]. 

The 5-dimensional gravity dual of a shock wave (ultrarelativistic nucleus) in our 4-dimensional 
space-time was first constructed in [6] . The AdSs shock wave metrics are shown below in Eqs. (|2.1| ) 



and ( |2.2j ). These metrics are solutions of Einstein equations in the AdSs bulk without sources. In 
four dimensions they correspond to nuclei of infinite transverse extent with a uniform distribution of 
matter in the transverse plane. Collisions of shock waves in Eqs. (12. If) and (|2.2j ) have been studied 
in [7,9,10,15] with the goal of explicitly constructing the metric after the collision. While the shock 
wave collisions in AdS3 allowed for an exact solution of the problem [7], it turned out to be signif- 
icantly harder in AdSs, allowing only for a perturbative solution of Einstein equations in graviton 
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exchanges [9,10]. While all-order graviton exchanges with one of the shock waves (corresponding 
to proton-nucleus collisions) were resummed exactly in [15], the full problem of nucleus-nucleus 
collisions involving all-order graviton exchanges with both nuclei still remains unsolved in AdSs. 

In [8] an alternative to the exact solution of Einstein equations was proposed: the authors of [8] 
constructed a trapped surface for a head-on collision of two shock waves with sources in the AdS 
bulk following [16, 17]. Sources in the bulk lead to the nuclei in the boundary theory having some 
transverse coordinate dependence in their matter distributions. Formation of a trapped surface 
before the collision indicates that a black hole will be formed in the future, after the collision. Thus 
the authors of [8] have proven that black hole is formed in a collision of two shock waves with point 
sources in the bulk. Generalizations of [8] to the case of nuclear collisions with non-zero impact 
parameter were presented in [11,12]. Also a trapped surface was found in [11] for an important case 
of collision of two shock waves with extended (not point-like) bulk sources. 

However, the exact implications of a source in the bulk for the boundary theory are still not 
entirely clear. The same energy-momentum tensor of the boundary theory can be given by metrics 
with extended sources at different bulk locations. It is possible that the sources would manifest 
themselves in fluctuations of the metric, but more research is needed to understand which bulk 
source gives the "right" fluctuations most accurately describing real-life heavy ion collisions. In [11] 
it was suggested that the position of the source in the bulk is related to the saturation scale of the 
shock wave. Initial steps on determination of saturation scale in shock waves were done in [18-20]. 
It appears more work is needed to clarify the complete impact of the bulk source on the boundary 
gauge theory. 

Interestingly the trapped surfaces found in [8,11,12] are always formed around the source in 
the bulk. One may therefore wonder whether the source is required for the trapped surface to form. 
No trapped surface analysis has been performed to date for the sourceless shock waves of Eqs. ( |2.1| ) 
and (|2.2| ) to answer this question. 

Here we perform a trapped surface analysis for a collision of two sourceless shock waves from 
Eqs. (|2.1| ) and (|2.2| ). We first consider the trapped surface obtained in [11] for a collision of two 
shock waves with extended sources in the bulk, and then take the limit in which the sources are 
moved to the deep infrared (IR) while keeping the energies of the shock waves in the boundary 
gauge theory fixed. Interestingly enough, the trapped surface does not disappear in this source-free 
limit, its lower boundary remains at finite value of the 5th dimension coordinate z with its finite 
area giving a finite expression for the produced entropy. We argue that collisions of two shock waves 
with sources in the deep IR (at z = oo) are indistinguishable from collisions of two shock waves 
without bulk sources by performing a perturbative solution of Einstein equations for the shock wave 
with sources in the bulk and taking the sources to z = oo. We also note that the trapped surface 
which remains after we send the sources to z = oo does not depend on how the limit was taken and 
on which sources were sent to infinity: the remaining trapped surface is the same for extended and 
point-like sources sent to the IR. 
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We therefore conclude that a collisions of two sourceless shock waves in AdSs leads to creation 
of a black hole in the bulk. The absence of bulk sources leaves no uncertainty in the interpretation of 
the physics and makes application of AdS/CFT correspondence better justified. For the boundary 
theory this result proves that thermalized quark-gluon plasma is produced in heavy ion collisions 
at strong coupling. 

The paper is structured as follows. In Sect. |2] we present the problem at hand and describe 
how the limit of sending the sources to the IR should be taken without changing the bulk physics. 
In Sect. |3] we present a lowest-order perturbative solution of Einstein equations for a collision of 
two shock waves with sources along the lines of a similar calculation for the sourceless shock waves 
in [10]. We take the limit of the shock waves sources going to the IR and show that our solution 
exactly maps onto the metric produced in a collision of two sourceless shock waves found in [10]. 
This provides a strong argument that the shock waves with sources at z = oo collide in the same 
way as the shocks without any sources. In Sect. [| we perform the trapped surface analysis and 
demonstrate that the trapped surface does not disappear when the sources are send to the deep IR. 
Thus we obtain the trapped surface for the collision of two sourceless shock waves. In Sect. [| we 
conclude by presenting a guess for the thermalization time inspired by our analysis (see also [9]). 
We argue that thermalization proper time is likely to be parametrically shorter than the light-cone 
stopping time for shock waves found in [10,15], which indicates that our conclusions may be applied 
to real-life heavy ion collisions at least at the qualitative level. We note however that the numbers 
generated by our approximate thermalization time estimate are too short to describe RHIC physics. 



2. The Problem 

High energy heavy ion collision can be realistically modeled by a collision of two ultrarelativistic 
shock waves. In [6], using the holographic correspondence [21], the geometry in AdSs dual to each 
one of the nuclei in the boundary theory is given by the following metric 

L 2 

ds 2 = — \-2dx + dx~ + t 1 (x~)z 4 dx' 2 + dx 2 ± + dz 2 } (2.1) 
z 2 ~~ 



for nucleus 1 and by 



L 2 

ds 2 = — {-2dx + dx' + t 2 (x + ) z 4 dx +2 + dx\ + dz 2 } (2.2) 

z 



for nucleus 2. Here dx\ = (dx 1 ) 2 + (dx 2 ) 2 is the transverse metric and x ± = (x° ± x 3 )/\/2 where 
x 3 is the collision axis. L is the radius of S5 and z is the coordinate describing the 5th dimension 
with the boundary of AdS 5 at z = 0. We have also defined 

ti(x-) = ^ (Zi_(iO), t 2 (x + ) = ^ (T 2++ (x + )} (2.3) 
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in accordance with the prescription of holographic renormalization [21]. Here (Ti {x~)) and 

(T 2++ (x + )) are the energy-momentum tensors of the two shock waves in the gauge theory. We 
assume that the nuclei are so large and homogeneous that one can neglect transverse coordinate 
dependence in (Ti (%~)) an d (T 2++ (x + )). Following [6] we take 

(Ti__(aT)> = M0O> (^++00) = MOO- (2.4) 

For simplicity we also put ^1 — ^2 — 

The metrics in Eqs. ( ^.1| ) and ( |2.2| ) solve Einstein equations in the empty AdSs space: 

1 6 

R P p - 2 9ti» R - 75 9i»> = °- ( 2 - 5 ) 

However, as we will see below, it is hard to perform the trapped surface analysis with the sourceless 
shock waves. To this end, as we have mentioned above, it will be more convenient to represent 
sourceless shock waves as limiting cases of the shock waves with sources, when the sources are sent 
to z = 00 while keeping energy-momentum tensor of the nuclei in the boundary theory intact. 

We therefore need to construct shock waves with sources in the bulk, which we will do following 
[8, 11]. We need to satisfy Einstein equations in AdS 5 with sources in the bulk 

1 6 

Rfw - g 9in> R - J2 9» u = % nG 5 Jfiv (2.6) 

where is the energy momentum tensor for bulk sources. We will not specify what fields contribute 
to create non-zero J^ u in the bulk: as for us the source will serve as an IR regulator we do not need 
to know the origin of in detail. The 5-dimensional Newton constant is 

G 5 = (2.7) 

C 

Eq. ( |2.6p can be rewritten as 

Rfiu + jp9nv = 8nG 5 ^J^ - J^j (2.8) 

with 

J = J/ = J^gT. (2.9) 

Following [11] for one shock wave we will consider a source without any transverse (a; 1 ,^ 2 ) 
coordinate dependence, with the only non-zero component of the energy-momentum tensor 

Ji°) = i-^-)^_ 2o ). (2.10) 
z L 
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The source is located at z = Zq and x~ = 0, spans the transverse directions and moves along the x + 
axis. E is a yet unspecified parameter with dimension of energy. To find the metric of the shock 
wave satisfying Eq. ( |2.8| ) with the source ( |2.10| ) we look for it in the following form generalizing 
Eq. O 

L 2 

ds 2 = — {-2dx + dx~ + <p(z)5(x')dx- 2 + dx 2 ± + dz 2 } . (2.11) 



Plugging Eqs. (|2.11|) and (|2.10| ) into Eq. ( |2.8|) we get the following equation for the " " component 



of Einstein equations [11] 



= 8 7rG 5 ^-S(z-z ). (2.12) 

A Z A Zq Li 



When solving this equation we require that <p(z) — > as z — > and that <p(z) is regular as z — > +oo. 
The latter condition is needed to avoid the singular behavior of metrics ( |2.1| ) and ( |2.2p in the IR. 
While the singularity of metrics ( |2.1| ) and ( |2.2| ) does not affect curvature invariants and is thus not 
unphysical, it is easier to perform trapped surface analysis which we intend to do below on a metric 
with is regular in the IR. 

Solving Eq. (|2.12|) with the boundary condition that <fi(z) — > as z — > and 4>(z) is regular at 
z — > +oo yields [11] 

;■--) = — r 1 - { (2-i3) 

1, Z > Z . 



Eqs. ( |2.13| ) and ( |2.11| ) give us the metric of a single shock wave with the bulk source ( |2.10| ) 



Using holographic renormalization [21] (see e.g. Eq. ( |2.3|) above) we conclude that the energy- 
momentum tensor corresponding to the metric ( |2.11| ) has only one non-zero component [11] 



(■/■__) = j^rS(x-) hm/M = ^<S (2.1 1) 

4 7TG5 2^0 2 4 Zq 



It is clear that this energy-momentum tensor would be the same as for the sourceless shock wave 
( ETTD given by Eq. ( ^4|) if we identify 

A* = ^ (2-15) 

obtaining 

(T__> = im5(x-). (2.16) 



The difference between the metrics for the shock wave with source in Eqs. ( |2.13| ) and ( p.ll|) 



and the sourceless shock wave in Eq. Q2.1|) is that the source regulates the metric in the IR. It is 
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important to note that if we take z — > oo limit of the metric in Eqs. fl2.ll ) and ( |2.13| ) keeping E/zq 
(and therefore u in Eq. ( |2.15|) fixed) we would recover the metric in Eq. ( |2.1|) without modifying 
the energy- momentum tensor of the gauge theory given by ( |2.16| ). At any finite z the metric of Eqs. 
( p. 11 ) and ( p. 13 ) becomes equivalent to (|2.1| ) in this limit, which sends the source at zq to the IR 
infinity. The question arises whether the metric ( |2TT| ) is equivalent to the zo —> oo, E/zq = const 
limit of the metric in Eqs. ( |2.11| ) and fl2.13| ). In other words, is having the sources at infinity 
identical to having no sources at all? 

We are interested in the answer to this question in the context of collisions of two shock waves. 
The question then becomes whether colliding shock waves from Eqs. (|2.1| ) and (|2.2|) are identical to 
colliding the shock wave in Eqs. ( |2.11|) and ( |2.13| ) with its counterpart with x + -H- x~ in the limit 
zo — > oo, E/zq = const of the resulting post-collision metric? 

The intuitive answer to the this question is "yes". Indeed it is highly unlikely that sources 
at zo = oo would affect any physics at finite z. Even in empty AdSs space light propagates with 
velocity 1 along the ^-direction. It would take light an infinite time to travel to any finite z from 
Zo = oo after the collision. The metric modification in the collision is only likely to lower the light 
velocity in the z-direction: in the "extreme" case when a black hole is created no signal from z = oo 
would be able to propagate outside of the horizon. Even more minor modifications of the metric 
are likely to only change the speed of light in ^-direction leaving it finite and not changing the 
above arguments. Hence any modification of sources at zo = oo in the collision is not going to affect 
the physics at finite z. Hence the collision of two shock waves with sources at zq = oo should be 
indistinguishable from the collision of two sourceless shock waves in Eqs. Q2.1J) and (|2.2|) . 

One may also think of a source at zo as providing an (externally imposed) infrared cutoff 1/zq 
on the transverse momenta kx of the partons inside the shock wave in the boundary gauge theory 
(see [11]). With this interpretation the limit of Zo —> oo, E/zq = const can be interpreted in the 
boundary theory as removing the IR cutoff on the transverse momenta of the partons while keeping 
the energy of the shock wave fixed. The shock waves without sources would then correspond to 
nuclei without an ad hoc IR cutoff on the transverse momenta of their partons in the boundary 
theory. Hence, from the standpoint of the boundary theory, the zo — > oo limit imposed on the 
four-dimensional shock waves dual to the shocks with sources in the bulk would simply remove the 
IR cutoff on partons' kj>. This would make the boundary theory shock waves identical to those 
dual to the sourceless shock waves in the bulk. Therefore, with the IR fc^-cutoff interpretation of 
1/zo [11] the zo — » oo limit also appears to be a justified way of obtaining duals of sourceless bulk 
shock waves in the boundary theory. 

To verify the above arguments we will perform a perturbative solution of Einstein equations 
for a collision of two shock waves with sources in the next Section. We will explicitly demonstrate 
that taking the Zq — > oo, Ejz\ = const limit of the obtained metric produced in the collision 
would simply reduce it to the metric produced in the collision of two sourceless shock waves found 
previously in [9,10], thus substantiating our intuitive argument above. 
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3. Perturbative Solution of Einstein Equations for Colliding Shock 
Waves with Bulk Sources 



Consider a collision of two shock waves with sources like the one given in Eq. fl2.10| ). The general 
metric for such a collision could be written as 

L 2 { 

ds 2 = — < — [2 + g(x + , x~ , z)] dx + dx~ + [(f>(z) 5(x~) + f(x + , x~~ , z)~\ dx~ 2 



+ 



[z)6(x + )+f(x + ,x-,z) 



dx 



+ 2 



[l + h(x + ,x ,z)] dx 2 ± + dz z 



(3.1) 



The functions /, /, g, and h are non-zero only for x + > 0, x~ > 0. Before the collision (for x~ < 
and x + < 0) the superposition of the metrics of colliding shocks (the terms with 0's above) solves 
Einstein equations Q2.8| ) exactly. 

We will follow [9,10,15] and find the functions /, /, g, and h perturbatively at the lowest order 
treating the shock waves as perturbations of the empty AdS 5 space. As <p(z) ~ fi one can argue 
that /, /, g, and h start at order fi 2 [10,15]. Our strategy is to expand Einstein equations to the 
order linear in /, /, g, and h and quadratic in 0. This is the same procedure as used in [10,15] for 
a collision of two sourceless shock waves. 

The main difference in the case at hand is that the shock waves now have sources. The energy- 
momentum tensors of the sources, given by the following non-vanishing components before the 
collision (order //, see Eqs. (|2.15| and ( |2.10|) )) 



J 



(o) 



/i jji 6(x )S(z-z ), J 



(0) 



Hj^5(x + )5(z-z ), 



(3.2) 



get modified in the collision. In principle to understand modifications of the bulk source one needs 
to know the field content of the source and the corresponding equations of motion for the fields. 
However, it turns out that this is not really necessary. Following a similar procedure for perturbative 
construction of classical Yang-Mills fields in nuclear collisions [22] we note that Einstein equations 
( UP imply 







(3.3) 



where is the covariant derivative. Imposing causality and using Eq. (|3.3| ) along with Einstein 
equations one can perturbatively construct the bulk energy-momentum tensor order-by-order in 
/i. Using the symmetries of the problem one can argue that it is unlikely that colliding sources 
would recoil in the transverse or z directions. This limits the non-zero contributions to the bulk 

energy-momentum tensor to J++, J and J + _ = J_ + . Note that to find at order fi 2 one only 

need the metric (|3.1|) at order \i. This means one does not yet need to know the functions /, /, g, 
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and h. It is then not too hard to infer the sources up to order /z 2 : the non- vanishing components 
of the bulk energy-momentum tensor are 

-v3 



J++ = p-±5{z- Zq) 

z 3 

J— = Vj^5(z- Zq) 



5(x + ) + i 6(x-) 5'(x+) [z <j)'{z) - <f>(z)] + . . . 
5(x~) + i 6(x + ) 5'(x-) [z (j)'(z) - 4>{z)\ + ... 



J-+ = —\l S(z - Zq) 5(x + ) 8(x ) 

L 6 



(z) + -z<f>'(z) 



+ ... 



(3.4a) 
(3.4b) 
(3.4c) 



Plugging Eqs. ( p. 1[ ) and ( |3.4| ) into fl2.8| ) and expanding the result in powers of \i we obtain at 
order /z 2 the following expressions for the "_L_L" and the u zz" components of Einstein equations 



(_U_) g z + 5 h z - z h zz + 2 z h x+ x - = 2 5{x + ) 5{x ) <j>(z) <j/(z) 

Zqh5(x + ) 6(x~) 5(z - z ) (f)'(z) 



16 7T G 5 ^ 5 
3 L 3 ~ 



(3.5a) 



(zz) g z + 2h z - zg zz - 2 z h z 



5(x + ) 5{x~) [-2 <j>(z) 4>'(z) + z {(t)'{z)f + 2z (f>{z) <f>"{z)] 



16 7T G5 

3 L 3 ~ 



— zlix8{x + )8(x )6{z - z )4>'(z). 



(3.5b) 



Here the subscripts indicate partial derivatives. Solving Eq. Q3.5a| ) for g z and substituting the result 
into Eq. ( |3.5b|) yields 



—3 h z + 3 z h z z — z 2 h z z z + 2 z 2 h x -\ 



167T G 



S(x + )5{x') 



5 „5 



3 L 3 



Eq. ( |3.6|) can be rewritten as 



^2 2 2 h x + x - 

z 



zlllZ [5'(Z - Z ) + 6(Z - Zq) cj>" \ Z )\ 



8(x+)8(x-) 



(3.6) 



x 



z[<P'(z)f 



16tt G 



^-Z^/IZ [5'(Z - Z ) (j)'(z) +5(z- Zq) <j>" {zj\ 



3 L 3 



(3.7) 



We can now substitute <p(z) from Eq. (|2.13| ) into Eq. Q3.7|) . There is a small subtlety: the derivative 
of <j>(z) is discontinuous at z — zq. It is therefore not clear which value of the derivative to choose, 
the one at z — zq — >■ + or the one at z — zq — > 0~. As for z > zq all derivatives of <fi(z) are zero, 
plugging the derivatives at z — z — > + into Eq. (|3.7| ) would simply eliminate all bulk source effects. 
It therefore seems more physical to use the derivatives at z — z — > 0~ . This gives 

2 



- h z — h zz + 2h x i 

z 



1 

3 



16 7T G 5 \i 

T 3 



5(x + ) 8(x 



-A8(z-z )-z 7 q6(z-Zq) 

(31 
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Eq. ( p.8| ) is easy to solve as the Green function for the operator on its left hand side was found 
in [15,23]. Defining the Green function by 



-d z -d 2 z + 2d + d- 



G(x + ,x ,z;x' + ,x' ,z') = S(x + — x' + ) 8(x — x' )5(z — z') (3.9) 
one can find an integral expression [15,23] 



G(x + ,x ,z;x' + ,x' ,z') = — 9(x + — x' + ) 9(x — x' ) — / dm 

L Z J 



x raJ (m a/2 (x + — x' + ) (x~ — x'~)^ J 2 (m z) J 2 (m z) 



which can be integrated to give 



G(x + ,x-,z-x' + ,x'-,z') = ^-9(x + -x' + )9(x- - x'~) 8(s) 8(2 - s) ,., 

2 71 * 1 



z l + 2s(s-2) 



(3.10) 



(3.11) 



with 



s = 



2 (x + - x' + ) (x~ -x'-)-{z- z') 

2~77' 



A2 



(3.12) 



With the help of Eq. Q3.10|) we solve Eq. (|3.8|) and write 

X + X oo oo 

h(x + , x~, z) = J dx' + J dx'~ J dz'-^—j J dram J Q (m\j2(x + — x' + )(x~ — x'~) j J 2 (m z) J 2 (m z') 

— oo — oo 

1 / 16 7T G 5 fi" 



L 3 



S(x' + ) 5(x'-) 



2 ■ 9(z 0(z' - zo) - z< 5(z' - z ) 



(3.13) 



Integrating over x' + and x' trivially yields 



oo oo 

z 2 



h(x + ,x~,z) = - I — — - — ] 9(x + )9(x~) I dz'-^— I dmm J (mr) J 2 (m z)J 2 (m z') 

3 V L 6 / I 2z' ' 



~/6 ~6 

Z — 9(z ~z')-f 9{z' - z ) - 4 5{z' - z ) 



where we defined the proper time 



V2 



X X 



(3.14) 



(3.15) 
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Let us evaluate the three terms in the brackets in Eq. ( |3.14|) separately. Start with the last 
term: it is proportional to 

J dz' ^—j J dmm J (mr) J 2 {m z)J 2 {m z') z 7 Q 5{z' — z ) = — Zq J dm m J (m t) J 2 (m z) J 2 (m z ) 
00 

1 + 2 s (s - 2) 



— z Q 9(s )9(2-so) 

2vr ^ So (2- So ) 



with 



s 



r 2 - (z-z Q )' 
2 z zo 



(3.16) 



(3.17) 



We see that taking zq 00 and keeping /1 fixed gives So ~ —z /(2z) such that the expression in 
Eq. (|3.16| ) becomes zero due to ^(so). Hence the last term in the brackets of Eq. (|3.14j ) does not 
contribute in the zq — >■ 00 limit. 

The second term in the brackets of Eq. ( p. 14 ) is proportional to 



dz'- 



z' 



dmm J (mr) J 2 {m z)J 2 (m z') 



2f) 



dm Jq {m r) J 2 (m z) J\ (m Zq) 







(3.18) 



2() 



with the last step being valid for zq > z + r, i.e., for the large zq we are interested in. 

We are left with the first term in the brackets of Eq. ( 3.14 ). Hence at large Zq we have 



h(x + ,x ,z) 



1 

3 

1 

3 



L 3 

8 7T G 5 \x 
L 3 



6(x + )6(x )z J dz' z j dmm J (mr) J 2 (m z)J 2 (m z') 
'0 



6(x + )6(x-)z 2 zt 



dm 



87iG 5 /u, 
L 3 



d(x + )d(x-)z 4 



2 1 2 

T H — z 
3 



J (mr) J 2 {mz) [6 J 4 (mz ) — rnz J 5 (mz )] 



(3.19) 



This is exactly the solution found for sourceless shock waves in [10]! Using h from Eq. (|3.19|) in 
Eq. ( |3.5a|) one would obtain function g, which, for z — )• 00 would also be Zo-independent and would 
also correspond to that found for sourceless shock waves in [10]. Similarly one can show that / 
and / would also reduce to the ones from [10] in the z — > 00 limit. We conclude that, at least at 
this lowest non-trivial order in /z, colliding shock waves with sources gives a metric which in the 
limit of zq — >■ 00 (keeping /i fixed) reduces to that produced in the collision of two shock waves 
without sources. This presents a strong argument supporting our earlier assertion that collisions of 
the shock waves with sources at zq = 00 are equivalent to collisions of the shock waves without the 
sources. 



- 10 - 



4. Trapped Surface Analysis 



Below we will present trapped surface analysis for a collision of two shock waves without bulk 
sources. We will begin by outlining general concepts of the trapped surface analysis and will 
present a naive attempt to find the trapped surface for a collision of shock waves from Eqs. (|2.1|) 



and ( |2.2| ). We will then obtain the trapped surface for a collision of two shock waves with bulk 
sources and take the limit of z — > oo, deriving the trapped surface for a collision of sourceless shock 
waves. We will solidify our above conclusion of the equivalence between the sourceless shock wave 
and the one with sources at z = oo by taking the limit of sources going to the IR for a collision of 
two different shock waves with extended sources at Z\ and z 2 and showing that the limiting trapped 
surface is the same as obtained before. 

4.1 Generalities 

Let us start with outlining some generalities of trapped surface. Consider the collision of two shock 
waves given by the following metric before the collision: 

L 2 L L 

ds 2 = — \-2dx + dx~ + dx\ + dz 2 } H $i(x±, z) S(x + ) dx +2 H ^(x±, z) S(x~) dx~ 2 (4.1) 

z l K J z z 



where (cf. Eq. (JOft) 

$i(x ± ,z) = -4>i(x ± ,z), z = l,2. (4.2) 



z 



The marginally trapped surface is found from the condition of vanishing of expansion 9 [24]. 
The trapped surface is made up of two pieces: S = S\ US2. S\{Si) is associated with shock wave at 
x + = (a; - = 0) before the collision. An additional condition is imposed requiring that the outer 
null normal to Si and S2 must be continuous at the intersection C = S\ D S2 point x + = x~ = to 
avoid delta function in the expansion. 

To calculate the trapped surface associated with shock wave at x + = 0, we use the following 
coordinate transformation [8,17]: 

x- -+x- + 9{x + ) (4.3) 

2 

to eliminate the delta-function discontinuity at x + = 0. 1 The trapped surface Si can then be 
parametrized by [17] 

x+ = 0,x~ = (4.4) 



1 Note a different definition for the light-cone coordinates used in [8,17]. 
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The condition of marginally trapped surface is the vanishing of expansion 9 = h^V with 
h^ u the induced metric and l v the outer null normal to the trapped surface. Similarly to [8,11,12], 
the condition gives rise to 

□ - -i J [^ 1 (x ± ,z)-^ 1 (x ± ,z)} = (4.5) 



with \l/i(a;_L, 2) = -ipi(x±, z) and the Laplacian is defined with respect to Euclidean AdS^ space 

ds 2 = % {dxl + dz 2 }. (4.6) 

z 1 

By analogy, we have the condition defining the trapped surface S 2 . 

(n - [* 2 (x ± ,z)-$ 2 (x ± ,z)} = 0. (4.7) 

The continuity of trapped surface Si and S 2 and their outer null normal on the cusp of the 
light-cone x + = x~ = reduce to the boundary conditions 

*i(x ± ,z)\c = Mx±,z)\c = (4.8a) 

V*i(x ±> z). W 2 (x ± ,z)| c = 8 (4.8b) 



where the boundary C is to be determined from Eq. (18). The covariant derivative V is again 
defined with respect to Eq. (|4.6|) . 

Having the equations for the trapped surface with arbitrary shock wave ( |4.5| ), fl4.7|) and ( |4.8| ) at 
hand, we are ready to apply them to the collision of source- free shock waves (|2.1| ) and ( j2.2| ). With 
the symmetry <f>i(z) = 4> 2 {z) = <p(z) (and thus ipi(z) = ip 2 (z) = ip(z)), they take a particularly 
simple form 

z 2 tf"(z) - z tf'(z) - 3 *(z) = (4.9a) 
*(Za) = *(z b ) = (4.9b) 

Z f 2 ^\z a f = ^'{z h f = 8. (4.9c) 

The boundary C in this case is given by z a < z < z&, as there is no dependence on transverse 
coordinates. Eq. ( |4.9|) is easily solved by 

tf( z ) = dz 3 + - (4.10) 

z 

with Ci and C2 arbitrary constants. 

Obviously we cannot have C\ = C 2 = because of Eq. ( |4.9c| ) . It is easy to see then Eq. ( [4.9b] ) 
would immediately require z a = zj,. Similar phenomenon of no trapped surface was observed in [17] 
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for collisions of gravitational shock waves in asymptotically Minkowskian 4- dimensional space-time. 
One may be tempted to conclude that trapped surface formation is not possible in collisions of 
source-free shock waves. Before accepting such conclusion, let us point out that the reason we 
choose C to be bounded by z a < z < z\> from both sides in the bulk is because the trapped surface 
has to be closed. However, AdSs is different from asymptotically Minkowskian spaces: it appears 
not quite clear whether the requirement of a closed trapped surface necessarily implies finite z b . 2 
If one searches for the trapped surface with z b = oo, i.e., with z > z a constraint only, such that 
conditions in Eqs. ( fOD are imposed only at z a , one gets 



giving 



L 

71 



4V2 



Z 



(4.11) 



(4.12) 



Unfortunately the conditions in Eqs. ( |4.9|) are insufficient to fix z a uniquely. 

However, z a in Eq. (|4.11|) can be fixed if we choose to study a closely relevant situation. Let us 
consider the trapped surface formation in the collision of two sourced shock waves 



ds 2 



ds 2 



T 2 

z 2 1 
T 2 

— \- 



-2dx + dx + 4>i(z) 5(x 
-2 dx + dx~ + 02 (z) 



dx 2 + dx\ + dz 2 } 
dx + 2 + dx\ + dz 2 } 



(4.13a) 
(4.13b) 



with the sources J 4 



^5{x 



e 2 l 2 



, - x .. , biz — zi) and J — , 

z\L V ' v l i Z2L 

the shock waves. As discussed in the previous sections, we keep ^Nr 
nuclei on the boundary have the same energy density. 

The equations (|4.9|) for the trapped surface now take the following form: 



— S(x ) 5(z — z 2 ) corresponding to each of 



z**»(z) 



-z%{z)-^i 
*i(z b ) = 



-16 7tG 5 Ei5(z — z^ 



fi such that the 



(4.14a) 
(4.14b) 



-^%(z b )%(z b ) 



(4.14c) 



where the boundary C is again z a < z < z b and i = 1, 2. Eq. ( [4.14| ) is solved by 



a 



D, 



■V/, 



Z < Zi 



Z > Zi 



(4.15) 



2 Requirement that the trapped surface has to be closed appears to stem from the cosmic censorship conjecture, 
which we assume to be true in AdSs x S§: the issue of whether trapped surfaces in AdSs necessarily have to be 
closed may require further investigation. 
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with the constants 



Ci 



D, 



4 3~~ 



(4.16) 



1 z a 



The third equation in ( 4.14Q gives the following simple relations: 

C\ C*2 



L 2 

Dl D 2 = -. 



4.2 Shock Waves with Identical Sources 



It is instructive to first consider a collision of identical shock waves in AdSs. Putting z\ 
and Ei = E 2 = E in Eqs. ( fl~T5l) and (|4~T6|) above we obtain from Eq. (E~P7D [11] 



Za + Z b = 



(4.17) 



z 2 — z 



(4.18) 



(Zq+Zj,) -3z Q Z b 
Z i Z 3 



1 

4- 



We want to take zq — >■ oo limit while keeping the energy of the shock wave in the boundary theory 
fixed. That is we want to hold 



EL 2 



(4.19) 



fixed. We rewrite Eq. ( |4.18 ) in terms of \i as 



z a + z b 



2 V2tt 



(Zg+Z b ) 2 -3 Za Zb 



(4.20) 



where we have replaced G§ = tt L 3 /2 iV 2 . Now, taking zq — » oo keeping /i fixed we can easily infer 
the asymptotics of z a and z\>. First one can consider the case that in this limit z a and z\> are of 
the same order, z a ~ zj,. In such case the first equation in Q4.20D gives z a ~ Zb ~ z q, which can 
not satisfy the second equation in (|4.2Q|) . As z a < Zb by definition, we are left to consider the case 
when, in the z — > oo limit one has z a <C z b . Then the first equation in (|4.20| ) yields 



2x/2 7r 2 

c 



n4 



(4.21) 
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which, when plugged into the second equation in ( |4.20|) along with the assumption that z a <C z b 



gives 



2a 



(4.22) 



The values of z a and Zb given by Eqs. (^4.22|) and ( |4.21| ) satisfy z a <C z b condition when z is 
large, which confirms that they give the correct asymptotics. In the strict Zq — > oo limit we see 
that Zb —> oo, but z a remains finite given by Eq. ( |4.22|) . Indeed Eqs. ( |4.22|) and (|4.21|) can also 
be obtained by solving Eqs. (fOC]) explicitly and taking the zo — > oo limit: the exact solution of 
Eq. ( |4.20[ ) giving real z a and z b is 



A 4 



x /2iV3e- 2W3 z^ + Az 8 p e 



r, ~4 



4? 



z b = + 77 V 211/3 £ 3 - 2 13 / 3 4 f + 4 2 8 /i 2 f 2 



(4.23a) 
(4.23b) 



where 



2V27 T 2 

iV 2 



(4.24) 



and 



,6 . /l 1 7 12 r.4 _ -12 .~2 



1/3 



(4.25) 



One can readily check that the 2:0 ~~ ► 00 asymptotics of Eqs. ( 4.23a ) and (|4.23b|) is given by Eqs. 
(p|) and (pi. 3 

Taking the zq — > 00 limit in Eq. ( |4.15|) one can see that the trapped surface is described by 



27T 2 



r 4 *4i 



V2 



[z 4 -fi-^] 



(4.26) 



which is exactly Eq. (|4.12|) with z a now fixed by Eq. ( |4.22| ). We see that introducing bulk source 
as a regulator of the metric in the IR and then taking zo — > 00 limit allows one to fix z a and hence 
determines the trapped surface uniquely. 

Now let us verify that the obtained value of z a in Eq. ( |4.22[) is independent of the way we take 
the limit of sending the bulk sources to infinite IR. Let us show that the same trapped surface arises 
in a more general case when the two shock waves are different from each other. 



3 Note that for zq < (2//2) 1 / 3 both z a and z\, from Eqs. (4.23a) and ( 4.23b| ) become complex and trapped surface 
ceases to exist: however this small- zq limit is the exact opposite of the zq — > 00 case we would like to consider here. 
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4.3 Shock Waves with Sources at Different Bulk Locations 



We consider a collision of shock waves with sources at different locations and with different ESs: 



now we have Z\ ^ z 2 and E\ ^ E 2 but with ^ L ^- 

z i 

A 2 and rewrite Eq. ( |4.17|) as 



-A 

z a z b 



E 2 L 2 

^2 



ji. It proves useful to set 



4 



Ai, 



I b 
72 T -J 

z b Z a 



_|_ l _ A!+A 2 + l 
A1A2 



.ft 

2(1 



2 r2 

z b a 



1 



27r 2 ^t y 2(l-AiA 2 ) 



(4.27) 



eliminating z x and z 2 - Finding solution for z a and z b seems to be a hard task. We instead first solve 
the first equation in ( |4.27| ) for z a /zb and use the obtained ratio in the second equation in ( |4.27| ) to 
find z a Zb. Using the product z a z b in = Ai, -5^ = A 2 we can write z\ and z% as (i = 1, 2) 



A, 



2 7T 2 /! 



4/3 



Ai + A 2 + l-3AiA 2 
2(l-AxA 2 ) 



2/3 [(A 1 + 1)(A 2 + 1)] 1 / 3 
Ai A 2 



(4.28) 



Again we have replaced G5 = n L /2 N 2 . We are interested in the limit Zx,z 2 — > 00 while keeping 
r\2 = 7* = finite and /x is fixed. It is not difficult to see that the limit can be achieved by taking 
Ai, A 2 — > 0. In this limit Eq. ( |4.28| ) takes a very simple form: 



1 


(f 


4A 2 




1 


(N 2 


4A a 





4/3 



4/3 



(4.29a) 
(4.29b) 



As z b > z a , the first equation in (|4.27| ) gives in the Ai, A 2 — > limit that z b ^> z a . Solving the second 
equation in ( (4.27[) for z b ^> z a one obtains z a asymptotics. Using the result in -^1 = Ai along with 
the first equation in (|4.27|) yields 



N 2 



1/3 



{ZlZ 2 ) 



—7- OO. 



(4.30a) 
(4.30b) 



These equations are completely analogous to Eqs. ([4.22|) and ( [4.2 1|) above. Therefore the trapped 
surface is independent of the way one send the bulk sources to the IR infinity: the sources do not 
have to be at the same bulk location to obtain the same answer as we had in the previous Subsection. 
To further test the independence of taking the limit of sources going to infinite IR bulk, we have 
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also taken a similar limit for the point-like sources, first advocated in [8] : the trapped surface found 
in [8] again reduced to the trapped surface found in this work above. 

This completes our analysis of trapped surface in the collision of two shock waves with sources 
infinitely deep in the bulk. We note unlike source in finite depth [11], no critical value for the 
energy density is found in the limit. The formation of the trapped surface is always guaranteed. 
The trapped surface is even independent of the ratio r%2 = i.e. the details of the limit! It is 
important to stress that the trapped surface does not disappear with the removal of the sources in 
the bulk, which can be viewed as IR regulators. In all the examples of scattering of shock waves 
with bulk sources the trapped surface always appears to be more or less centered around the source 
in the x±, z space. One was tempted to conjecture therefore that the trapped surface is an inherent 
property of non-zero bulk energy-momentum tensor. Our result proves otherwise, giving an example 
of the source-free shock waves collision with a well defined trapped surface. 

4.4 Limiting Trapped Surface 

To summarize our trapped surface analysis let us re-state that the profiles of the trapped surface 
are given by 

2tt 2 L/i , 3 ( z 3 z\ 



* M = ^* (431) 

which is exactly Eq. (|4.11|) with z* a from Eq. fl4.22|) . In the transformed light cone coordinates (see 
Eq. Q4.3Q) the trapped surface is then determined by 



*+ = 0,," = -^ [z* -***] = [z* - /T 4/3 ] (4.32) 



with an analogous expression for the other shock wave obtained by interchanging x + -H- x in 
Eq. (p|). 

The trapped surface for a collision of source-free shock waves from Eq. (|4.32 ) is illustrated in 



Fig. [I]. One can clearly see that the trapped surface is present at all times before the collision and 
rises from the deep IR toward finite values of z. Similar behavior was observed for the trapped 
surface in the numerical model of heavy ion collision involving gravitational perturbations in the 
4-dimensional world in [25,26]. A horizon rising from the deep IR was also deduced in [27] for a 
model of heavy ion collision involving a rapidity-independent matter distribution after the collision. 

It is interesting to point out that the obtained shape of the trapped surface appears to imply 
that the black hole produced in the collision would have a singularity at z = oo with the horizon 
independent of the transverse coordinates x±. This is indeed very similar to the black hole dual to 
Bjorken hydrodynamics constructed in [6]. The main difference is that in our case the metric (and 
the energy- momentum tensor in the gauge theory) are rapidity-dependent, as follows from explicit 
calculations of the metric produced in shock wave collisions [9, 10, 15]. 
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shock wave 1 collision shock wave 2 
position 



Figure 1: An illustration of the trapped surface in the collision of two sourceless shock waves. Vertical 
axis is the bulk z-direction, the horizontal left-right axis can be thought of either as the collision axis or 
as the time direction. The trapped surface is shaded. 



Our estimate for the produced entropy per unit transverse area A± for a collision of two shock 
waves with the sources at z is [11] 



S N% 



A] 



2tt 



zi zi 



(4.33) 



Using Eqs. fl4.22Q and ( |4.21| ) we obtain for z — > oo 



S 



2 ..21 1/3 



As /i 2 ~ s with s the center of mass energy of the collision, we get 



(4.34) 



_S_ 
XT 



oc s 



1/3 



(4.35) 



in agreement with the result obtained in [8]. 

The entropy from Eq. ( |4.33[ ) is plotted in Fig. [2] in arbitrary units as a function of the bulk 
source location z (for a collision of two identical shock waves). Fig. |2] demonstrated that produced 
entropy becomes practically independent of the bulk source position rather fast, approaching its 
asymptotic value well before zq jl 1 ^ 3 becomes large. 

As we noted above, for zq jl 1 ^ 3 < 2 1 / 3 both z a and z& given by Eqs. fl4T23aD and (|4.23b|) become 
complex and the trapped surface ceases to exist (see also [13] for a similar result). This likely 
implies that no black hole is formed in collisions of such shock waves. To understand this result 
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1.5 2.0 2.5 3.0 

Figure 2: The (lower bound on the) entropy density produced in the collision of two identical shock waves 
with sources as a function of the source position z$ in units of /2 -1 / 3 . The entropy density is in arbitrary 
units. 



from the boundary gauge theory perspective one has to have a rigorous interpretation of what 
shock wave sources in the bulk are dual to in the gauge theory. Such interpretation is missing at 
the moment, which inspired our present investigation of collisions of the sourceless shock waves. 
We may speculate though: following [11] we may assume that the inverse position of the source in 
the bulk 1/zo provides an IR cutoff on the transverse momenta of the partons in the shock waves' 
wave functions in the boundary theory. Reducing zq would increase the cutoff 1/zq thus decreasing 
the number of partons: this is likely to lower the number of degrees of freedom produced in the 
collision, leading to the reduction of the entropy density with decreasing zq in Fig. ^|. Still it is not 
entirely clear why the trapped surface disappears completely at a finite small zq forcing the estimate 
for produced entropy to go to zero. Indeed our delta-function shock waves are described by a single 
dimensionful parameter fx (or /i): the largest momentum scale in the problem is therefore /i 1//3 . If 
1/zq is the IR cutoff, then clearly it can not exceed the largest momentum scale: hence 1/zq < /i 1//3 . 
This, however, can not explain why the trapped surface vanishes entirely at z = 2 1 / 3 /t -1 / 3 . Besides 
nothing pathological seems to happen in the perturbative solution presented in Sect. |3| for small 
finite zq. In this work we are interested in the large- zq asymptotics: Fig. |2] demonstrates that the 
produced entropy density does not seem to change much between having sources at finite large 
Zq > 2 1 / 3 /} -1 / 3 and having no sources at all, which seems to agree with the IR cutoff interpretation 
of the sources and, more importantly, shows that the entropy is "well-behaved" in the Zq — > oo limit 
we are taking. We leave the detailed study of the small- zq regime for future work. 
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5. Thermalization Time Estimate and Conclusions 



The result fixing z* in Eq. (|4.22|) could be predicted if one realizes that in the limit of delta-function 
shock waves the problem has only one dimensionful parameter p, which has dimensions of mass 
cubed. If a non-vanishing trapped surface is created in such collisions it has to be proportional to 
the only distance scale in the problem: 1/fl 1 ^ 3 . Stretching this analogy further one should expect 
that the proper time of thermalization (the time of black hole formation) is 

(5.1) 

as was originally suggested in [9]. 

An interesting question is the relation between this thermalization time and the time it takes 
for the shock waves to stop. It was argued in [10,15] that colliding shock waves come to a complete 
stop shortly after the collision. One can argue that fi ~ p + A 2 A 1 ^ 3 [10], where p + is the large 
longitudinal momentum of a "nucleon" in the shock wave, A is the typical transverse momentum 
scale in the shock, and A is the atomic number of the nucleus we model by the shock wave. The 
characteristic light-cone stopping time for a shock wave moving in the light-cone "plus" direction 
is given by [10, 15] 

,1 . . 

x sto P ~ (5-2) 

This is of course parametrically much longer than 

11 

Ttherm ~ ^ ~ ( p + A 2 A 1 / 3 ) 1 / 3 ' ^ 

Hence, if one assumes that thermalization happens at mid-rapidity first, then, as near mid-rapidity 
t ~ r, the time of thermalization is Uherm ~ Ttherm "C t stop = x^ top /\/2. It is therefore likely that 
thermalization happens at times which are parametrically earlier than the stopping time. If our 
guess of thermalization time is correct, this would imply that the shock waves still move along their 
light cones when thermalization happens, justifying an assumption commonly used in hydrodynamic 
simulations of heavy ion collisions. Note also that the thermalization time in Eq. ( |5.3|) is very 
short, and decreases with center-of-mass energy of the collision. (In fact, as was noticed in [9] this 
thermalization time is too short: if one plugs in p + = 100 GeV, A = 0.2 GeV and A = 196 into 
the parametric estimate ( |5.3| ) one would obtain Ttherm ~ 0.07 fm/c for RHIC, which is far too short 
for agreement with hydrodynamic simulations [1,2]. Indeed the thermalization time estimate of 
Eq. ( p.lfl is too crude for 0.07 fm/c to be taken literally, and a numerical coefficient in front of 
the estimate ( |5.1|) , if it results from a more exact calculation and from a more realistic model of 
colliding nuclei, may significantly change this number.) 

It is important to stress the difference between the mathematical limit of delta-function shock 
waves (a — > with a being the x~-width of the smeared non-delta-function shock wave [10, 15] 
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moving in the x + direction or vice versa) and the physical high energy limit of p + ^> A for nuclei. 
While in the former limit Jx is the only non-vanishing dimensionful parameter in the problem, the 
latter limit has another non- vanishing dimensionful scale jx a, which in fact gives the stopping time 
(f5.2|), x~^ top ~ [10,15]. (As one can easily see a ~ A 1 ' 3 /p + in the center-of-mass frame, such 

that jla ~ A 2 A 2 ! 3 is independent of p + [10, 15].) With the presence of two momentum scales in 
the problem the validity of the thermalization time estimate of [9] shown here in Eq. ( |5.1[ ) becomes 
less apparent. Our trapped surface analysis resulting in Eq. (|4.22|) appears to indicate that it is 
the momentum scale which depend only on ft and not on a that matters for thermalization, thus 



providing new evidence to support the estimate in Eq. (|5.1| ). In other words we show that if one 
neglects the smaller second momentum scale ft, a and approximates the shock wave profiles by delta- 
functions, thermalization is achieved in the collisions at the time given in Eq. ( |5.1|) . If one treats the 
problem more carefully and includes the scale fl a by considering shock waves of finite longitudinal 
spread [10,15], Eq. ( f5.1|) is likely to get corrections with the relative suppression factor being some 
positive power of jla/fl 2 ^ 3 , which is very small for high energy collisions, thus leaving the estimate 
in Eq. Q5.1] ) practically unchanged. 

One may argue that the strongly-coupled dynamics of the J\f = 4 SYM medium produced in 
shock wave collisions may be similar to that of strongly-coupled QCD medium. Then our conclusion 
of rapid thermalization may be applicable to soft (non-perturbative, kr ~ ^qcd) modes in heavy ion 
collisions, which would thermalize very quickly. Harder (perturbative) modes may then thermalize 
through interactions with the soft non-perturbative thermal bath, though more work is needed to 
justify such thermalization scenario and to modify the thermalization time estimate ( [5.3D to take 
into account perturbative dynamics. 

Another interesting question would concern understanding the relation between the rather 
quick thermalization in heavy ion collisions for the theory at strong coupling argued here and the 
impossibility of thermalization at weak coupling suggested in [28] by one of the authors. 4 While 
further research is needed to clarify this problem, the solution may have already been suggested 
in [32,33], where the authors argue that it is possible that there is a critical value A c of 't Hooft 
coupling A. For A > A c black hole formation is likely in high energy collisions. At the same time, 
for A < A c the black hole is not formed in high energy collisions [32,33], corresponding to no 
thermalization in the boundary theory. Indeed in the case of real-life heavy ion collisions, due to 
the running of the strong coupling constant, the coupling assumes a wide range of values in a single 
collision. The coupling is always large for soft transverse modes, making thermalization due to large 
coupling effects likely in light of our above results. 

To conclude let us point out once more that we have obtained a trapped surface for a collision 
of two sourceless shock waves in AdS 5 . The shape of the trapped surface is given by Eq. (|4.32| ) and 
is illustrated in Fig. [I]. Existence of this trapped surface proves that a black hole is created in the 
bulk for a collision of two sourceless shock waves, corresponding to creation of thermalized medium 



^ote that perturbative thermalization scenarios have been advocated in [29-31]. 
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(quark-gluon plasma) in the boundary gauge theory. 
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